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^ ■ Abstract 

a. 

This paper focuses on the equidimensional decomposition of affine vari- 
eties defined by sparse polynomial systems. For generic systems with fixed 
supports, we give combinatorial conditions for the existence of positive di- 
mensional components which characterize the equidimensional decomposition 
of the associated affine variety. This result is applied to design an equidi- 
mensional decomposition algorithm for generic sparse systems. For arbitrary 
^r) ■ sparse systems of n polynomials in n variables with fixed supports, we obtain 

an upper bound for the degree of the affine variety defined and we present 
an algorithm which computes finite sets of points representing its equidimen- 
sional components. 



$ '• 1 Introduction 



The aim of this paper is to describe the affine solution set of a polynomial system 
taking into account the sets of exponents of the monomials with nonzero coefficients 
in the polynomials involved, that is, their support sets. 

Bernstein [T], Kushnirenko [2T] and Khovanskii [TH] proved that the number of 
isolated solutions in (C*) n of a polynomial system with n equations in n variables 
is bounded by a combinatorial invariant (the mixed volume) associated with their 
supports. This result, which may be considered the basis for the current study 
of sparse polynomial systems, hints at the fact that the algorithms solving these 
systems should have shorter computing time than the general ones. 

There are several algorithms to compute either numerically or symbolically the 
isolated roots of sparse polynomial systems in (C*) n (see, for example, [33], [i~4"] . 
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[26J, [ 1 TJ ) . The efficiency of some of these algorithms relies on the use of polyhedral 
deformations preserving the monomial structure of the polynomial system under 
consideration. 

The first step towards the study of the solutions of sparse systems in the affine 
case was to obtain upper bounds for the number of isolated solutions in C n in terms 
of the structure of their supports and to design numerical algorithms to compute 
them (see [25], [23], [27], [15], [6], [TJ). Symbolic algorithms performing this task 
were given in [T7] and [T3] . 

The following natural step is to characterize the components of higher dimension 
of the affine variety defined by a system of sparse polynomial equations taking into 
account their supports. In this context, in |32j, certificates for the existence of curves 
are given in the numerical framework. 

There are different symbolic algorithms describing the equidimensional decompo- 
sition of a variety which only take into consideration the degrees of the polynomials 
defining it and not their particular monomial structure. The most efficient deter- 
ministic ones can be found in [5] and [5]. Probabilistic algorithms with shorter 
running time are given in [5] and [18]. The complexities of these probabilistic algo- 
rithms are polynomial in the Bezout number of the system, which, in the generic 
case, coincides with the degree of the variety the system defines. Other probabilistic 
algorithms are presented in [22] and [IB] with complexities depending on a new in- 
variant related to the system (the geometric degree) which refines the Bezout bound. 
From the numerical point of view, equidimensional decomposition algorithms based 
on homotopy continuation are proposed in [30] . 

In this paper we analyze, both from the theoretic and algorithmic points of 
view, the equidimensional decomposition of the affine variety defined by a sparse 
polynomial system. 

First, we consider the case of generic sparse systems. In this context, there ex- 
ists a major difference with the case of dense polynomials. The set of solutions of 
a generic system of n polynomials in n variables with fixed degrees consists only of 
isolated points. However, fixing the set of supports of the n polynomials in n vari- 
ables involved in a sparse system, for generic choices of its coefficients, there may 
appear affine components of positive dimension (see, for instance, Examples [T] and [2] 
below). We show that the existence of these generic components of positive dimen- 
sion depends only on the combinatorial structure of the supports: in Proposition 
H] below, we give conditions that yield these components. Such conditions provide 
not only a theoretic description of the equidimensional decomposition of the affine 
variety V(f) defined by a generic sparse system f in terms of the solution sets in the 
torus of smaller systems fj associated to subsets / C {1, . . . ,n} but also a formula 
for the degree of ^(f) in this generic case (see Theorem [5] below) . 

This result is used to design a probabilistic algorithm which, for a generic sparse 
system f, computes the equidimensional decomposition of V(f) with a complexity 
depending on its degree and combinatorial invariants associated with the system 
(see Theorem [8]). The idea of the algorithm is to compute first a family of subsets 
/ C {1, . . . , n} which may lead to components of V(f) and solve the corresponding 
polynomial systems fj by applying symbolic polyhedral deformations ( [1 TJ ) and a 
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Newton- Hensel based procedure QllJ). The output of the algorithm is, for each 
k = 0, ... ,n — 1, a list of geometric resolutions representing the equidimensional 
component of dimension k of V(f). A geometric resolution of an equidimensional 
variety is a parametric type description of the variety which is widely used in sym- 
bolic computations (see, for instance, [TO], [HJ, i n Section I2TTI below we 
give the precise definition we use. 

The following step is to consider the equidimensional decomposition of the affine 
variety defined by an arbitrary system of sparse polynomials. A question to answer 
beforehand is which parameter should be involved in the algebraic complexity of 
an algorithm solving this task. From previous experience, a natural invariant ex- 
pected to appear in the complexity bounds is the degree of the variety, which is, in 
particular, an upper bound for the number of its irreducible components. 

Unlike the Bezout bound for dense polynomials, in the sparse setting, the degree 
of the affine variety defined by a generic square system is not an upper bound for the 
degree of the variety defined by any system with the same supports (see Example 
E]). In |20j, a bound for the degree of the affine variety defined by an arbitrary 
sparse polynomial system depending on a mixed volume related to the union of the 
supports of the polynomials is presented. Here, we obtain a sharper bound for this 
degree also given by a mixed volume associated to the supports but not involving 
their union: 

Theorem. Let f = (fi,...,f n ) be n polynomials in C[Xi,...,X n ] supported on 
A = (A 1: ...,An) and let V(f) = {iG C n | fi(x) = for every 1 < i < n}. Then 

deg(V(f)) < MV n (Ai U A, . . . , A n U A), 

where A = {0, e±, . . . , e n } with Ci the ith vector of the canonical basis of W 1 and 
MV n stands for the n-dimensional mixed volume. 

Finally, we obtain an algorithm which, using a polyhedral deformation, describes 
points in every irreducible component of the affine variety defined by an arbitrary 
square sparse system with complexity depending on the degree bound previously 
stated. 

The idea of the algorithm relies on the fact that cutting the variety with a generic 
affine linear variety of codimension k, sufficiently many points in each irreducible 
component of dimension k can be obtained. To keep the complexity within the 
desired bounds, instead of computing this intersection, we proceed in a particular 
way which enables us to compute a finite superset of the intersection included in 
the variety. Representing a positive dimensional variety by means of a finite set of 
points is a well-known approach in numerical algebraic geometry (see the notion of 
witness point supersets in |30j). 

Theorem. Let f = (f\, . . . , f n ) be n polynomials in Q[Xi, . . . , X n ] supported on A = 
(Ai, . . . , An). There is a probabilistic algorithm which, taking as input the sparse 
representation of f computes a family ofn geometric resolutions (R^°\ . . . , R^ n ~^) 
such that, for every < k < n — 1, R^ k ' represents a finite set containing deg Vfc(f) 
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points in the equidimensional component 14 (f) of dimension k ofV(f). The com- 
plexity of the algorithm is of order 0~(n 4 dND 2 ) , where d = maxi<j< n {deg(/j)}, 
N = £" =1 #(Aj U A) and D = MV n (Ai U A, . . . , An U A). 

Here 0~ refers to the standard soft-oh notation which does not take into account 
logarithmic factors. Furthermore, we have ignored factors depending polynomially 
on the size of certain combinatorial objects associated to the polyhedral deformation. 
For a precise complexity statement, see Theorem [TT1 below. 

The paper is organized as follows. In Section the basic definitions and nota- 
tions used throughout the paper are introduced. Section [3] is devoted to the equidi- 
mensional decomposition of affine varieties defined by generic sparse systems: first, 
we consider the solution sets in (C*) n of underdetermined systems (see Subsection 
13. ID ; then, we prove our main theoretic result on equidimensional decomposition 
and present our algorithm to compute it (see Subsection 13.21) . Finally, in Section 
HJ we consider the case of arbitrary sparse systems: we prove the upper bound for 
the degree of affine varieties defined by these systems in Subsection 14.11 and, in 
Subsection I4.2[ we describe our algorithm to compute representative points of the 
equidimensional components. 

2 Preliminaries 

2.1 Basic definitions and notation 

Throughout this paper, we deal with polynomials in Q[X l5 . . . ,X n ], that is to say 
polynomials with rational coefficients in n variables X = (Xi, . . . , X n ). If f = 
(/i> • • • j fs) is a family of such polynomials, V(f) will denote the algebraic variety 
of their common zeroes in C ra , the n-dimensional affine space over the complex 
numbers. 

The algebraic variety V(f) C C n can be decomposed uniquely as a finite union 
of irreducible varieties in a non redundant way. This leads to the equidimensional 
decomposition of the variety: 

n 

v{f) = \Jv h (f), 

k=0 

where, for every < k < n, 14(f) is the (possibly empty) union of all the irreducible 
components of dimension k of V(f). 

The degree of each equidimensional component 14(f) is the number of points in 
its intersection with a generic affine linear variety of codimension k, and the degree of 
V(f), which we denote by deg(V(f)), is the sum of the degrees of its equidimensional 
components (see [T2]). 

A common way to describe zero-dimensional affine varieties defined by polyno- 
mials over Q is a geometric resolution (see, for instance, [TTJ and the references 
therein). The precise definition we are going to use is the following: 
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Let V = . . . ,^ D ^} C C n be a zero-dimensional variety defined by rational 
polynomials. Given a linear form £ = £\X\ + • • • + £ n X n in Q[Xi, . . . , X n ] such that 
^ £(^) if i 7^ j, the following polynomials completely characterize V: 

• the minimal polynomial q = Yli<i<D( u ~ )) Q[ u ) °f @ over t ne variety 

V (where u is a new variable), 

• polynomials v\, . . . , v n G Q[w] with deg(fj) < -D for every 1 < j < n satisfying 

V = {(vifa), . . . , 1^(77)) G C n I 77 G C, g(r/) = 0}. 

The family of univariate polynomials (q,v\, . . . ,v n ) G Q[w] n+1 is called a geometric 
resolution of V (associated with the linear form £). 

An equivalent description of V can be given through the so-called Kronecker 
representation (see [H]), which consists of the minimal polynomial q and polyno- 
mials wi,...,w n e Q[u] such that V = {(^(77), • • • , ^(v)) G C n | 77 G C, 4(77) = 0}. 
Either representation can be obtained from the other one in polynomial time. 

The notion of geometric resolution can be extended to any equidimensional va- 
riety: 

Let V C C™ be an equidimensional variety of dimension r defined by poly- 
nomials fi, . . . , f n - r G Q[X]. Assume that for each irreducible component C of 
V, the identity 1(C) D Q[Xi,...,X r ] = {0} holds. Let i be a nonzero linear 
form in Q[X r+1 , . . . , X n ] and 1x1 : V — > C r+1 the morphism defined by irg(x) = 
(xi, . . . ,x r ,£(x)). Then, there exists a unique (up to scaling by nonzero elements 
of Q) polynomial Qt G Q[A 1; . . . , X r , u] of minimal degree defining 7ii(V). Let 
qt G Q(Xi, . . . , X r )[u] denote the (unique) monic multiple of Qi with deg^g^) = 
deg u (Qe). In these terms, if £ is a generic linear form, a geometric resolution of V 
is (qe, v r+ i, . . . , v n ) G (Q(Xi, . . . , X r )[u)) n ~ r+1 , where, for r + 1 < i < n, Vi satisfies 

^(£) X t = v t {£) in Q(X U X r ) ® Q[V] 
and deg M (?j i ) < deg u (q e ). 

2.2 Algorithms and codification 

Although we work with polynomials, our algorithms only deal with elements in Q. 
The notion of complexity of an algorithm we consider is the number of operations 
and comparisons in Q it has to perform. We will encode multivariate polynomials 
in different ways: 

• in sparse form, that is, by means of the list of pairs (a, c a ) where a runs over the 
set of exponents of the monomials appearing in the polynomial with nonzero 
coefficients and c a is the corresponding coefficient, 

• in the standard dense form, which encodes a polynomial as the vector of its 
coefficients including zeroes (we use this encoding only for univariate polyno- 
mials), 
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• in the straight-line program (sip for short) encoding. A straight-line program is 
an algorithm without branchings which allows the evaluation of the polynomial 
at a generic value (for a precise definition and properties of sip's, see [2]). 

In our complexity estimates, we will use the usual O notation: for /, g : Z> — > K, 
f(d) = 0(g(d)) if | /(d) | < c|g(<f)| for a positive constant c. We will also use 
the notation M(d) = d\og 2 (d) log(log(c?)), where log denotes logarithm to base 2. 
We recall that multipoint evaluation and interpolation of univariate polynomials of 
degree d with coefficients in a characteristic-0 commutative ring R can be performed 
with 0(M(d)) operations and that multiplication and division with remainder of 
such polynomials can be done with 0(M(d) / log(d)) arithmetic operations in R. 

We denote by Q the exponent in the complexity estimate 0(d Q ) for the multipli- 
cation of two d x d matrices with rational coefficients. It is known that Q < 2.376 
(see [SI Chapter 12]). Finally, we write Q for the exponent (Cl < 4) in the complexity 
0(d n ) of operations on d x d matrices with entries in a commutative ring R. 

Our algorithms are probabilistic in the sense that they make random choices of 
points which lead to a correct computation provided the points lie outside certain 
proper Zariski closed sets of suitable affine spaces. 

2.3 Sparse systems 

Given a family A = (Ai, ■ ■ ■ ,A S ) of finite subsets of (Z> ) n , a sparse polynomial 
system supported on A is given by polynomials fj = ^2 aeAj Cj ta X a , in the variables 
X = (X\, . . . ,X n ) with Cj a 6 C \ {0} for each a G Aj and 1 < j < s. We write 
f = (ft, . . . , f s ) for this system. 

Assume s = n. We denote by MV n (A) the mixed volume of the convex hulls of 
A±, . . . ,A n in M. n (see, for example, [U Chapter 7] for the definition), which is an 
upper bound for the number of isolated roots in (C*) n of a sparse system supported 
on A (see [T]). 

The mixed volume MV n (A) can be computed as the sum of the n- dimensional 
volumes of the convex hulls of all the mixed cells in a fine mixed subdivision of A. 
Such a subdivision can be obtained by means of a standard lifting process (see [HI 
Section 2]): let lo = (ui,...,uj n ) be a n-tuple of generic functions Uj : Aj — >■ M 
and consider the polytope P in IR n+1 obtained by taking the pointwise sum of the 
convex hulls of the graphs of Uj for 1 < j < n. Then, the projection of the lower 
facets of P (that is, the n-dimensional faces with inner normal vector with a positive 
last coordinate) induces a fine mixed subdivision of A. The dynamic enumeration 
procedure described in [21] appears to be the fastest algorithm known up until now 
to achieve this computation of mixed cells. 

The stable mixed volume of A, which is denoted by SM n (A), is introduced in 
[j~5] as an upper bound for the number of isolated roots in C n of a sparse system 
supported on A. 



6 



3 Generic sparse systems 



3.1 Toric components 

Let n and m be positive integers and let A = (A\, . . . ,A m ) be a family of finite 
subsets of (Z> ) n . For 1 < j < m, let 

Fj(Cj,X) = ^ C iA Xa 

where X = (Xx,...,X n ); for a = (ax,...,a n ), X a = Yli<j< n ^ 3 , and Cj = 
(Cj,a)aeAj are Nj = i^Aj indeterminate coefficients. 
Following [31], consider the incidence variety 

{(x, c) G (C*) n x (P^- 1 x ■ • • x P^-- 1 ) | Fj(cj, x) = for every 1 < j < m} 
and its projection to the second factor 

Z = {ce P^^x- ■ •xP w '^ 1 | 3 a; G (C*) n with Fj( Cj ,x) = for every 1 < j < m). 

Note that the elements in Z correspond to the coefficients of the systems sup- 
ported on A which have a solution in (C*) n . 

Lemma 1 The Zariski closure of Z equals p-^ 1-1 x ■ • • x P Wm_1 if and only if, 
for every J C {l,...,m} ; dim (j2j eJ Aj j > #J. In particular, if m > n, a 
generic system supported on A has no solutions in (C*) n . Moreover, if m < n 
and dim Aj^j > #J for every J C {1, . . . ,m}, the solution set in (C*) n of a 

generic system supported on A is an equidimensional variety of dimension n — m and 
degree MV n (Ai, • • • , Am, A^ n_m ^) ; where A = {0, e±, . . . , e n } with e« the ith vector 
of the canonical basis of IR n and the superscript (n — m) indicates that it is repeated 
n — m times. 

Proof: The first statement of the Lemma follows as in [3T| Theorem 1.1]. 

Assume that m < n and for every J C {1, . . . , m}, dim (^2j e J Aj^j > # J. Then, 
if A = (A, • • • , A m , A( n ~ m )), we have that for every J C {1, . . . , n}, the inequality 
dim Aj^ > #J holds, since dim(A) = n. Therefore, MV n (A) > 0, which 

implies that a generic system supported on A has finitely many solutions in (C*) n 
(as many as MV n (A)). 

Now, the solution set of a generic system supported on A is the intersection 
of the solution set of a generic system of m equations supported on A, which is a 
variety of dimension at least n — m in (C*) n , and n — m generic hyperplanes. We 
conclude that the solution set in (C*) n of a generic system supported on A is an 
equidimensional variety of dimension n — m and degree MV n {A) . □ 
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Assume now that f = (fx, ■ ■ ■ , fm) are generic polynomials in the variables 
X = (X h ...,X n ) supported on A = (Ai, Am) C (Z^ ) m , with m < n. The 
previous lemma states that the affine variety V*(f) C C n consisting of the union of 
all the irreducible components of V(f) = {x G C n | fj(x) = for every 1 < j < s} 
that have a non-empty intersection with (C*) n is either the empty set or an equidi- 
mensional variety of dimension n — m. 

In what follows, we extend the symbolic algorithm from [T71 Section 5], which 
deals with the case m = n, to a procedure for the computation of a geometric 
resolution of V*(f) for arbitrary m < n. As in [17], our algorithm assumes that 
a fine mixed subdivision of (A, A^ n ~ m ^) induced by a generic lifting function u = 
(cji, . . . ,cj n ) is given by a pre-processing whose cost will not be considered in our 
complexity estimates. 



Algorithm GenericToricSolve 

INPUT: A sparse representation of a generic system f = (fx, ... , f m ) in the vari- 
ables X = (Jfi, . . . , X n ) supported on A = (Ai, . . . , A m ), a lifting function uj = 
(ui, . . . , uj n ) and the mixed cells in the subdivision of (A, A( n ~ m ^) induced by to. 

1. If the fine mixed subdivision of (^4., A^" - "^) does not contain any mixed cell, 
return R = 0. Otherwise, continue to Step [2j 

2. Choose randomly the entries of a matrix A = (clm) G Q nxn and a vector 
6= (&!,...,&„) GQ". 

3. For 1 < h < n — m, consider the affine linear forms = Ym=i a hi^i ~ 

4. Apply [TT|, Algorithm 5.1] to obtain a geometric resolution (q(u), Vi(u), . . . , v n (u)) 
of the isolated common zeroes of f , L±, . . . , L n _ m in (C*) n . 

5. Obtain an sip for the polynomials g := f(A~ 1 Y) in the new variables Y = 
(Y h ...,Y n ). 

6. Compute (wx(u), . . . , w n (u)) t := A(vi(u), . . . , u n («))*. 

7. Apply the Global Newton Iterator from [TTJ Algorithm 1] to the polynomials 
g(Y), the geometric resolution (q(u),w n - m+ i(u), . . .,w n (u)) of V(g(b, Y n -m+i, 
and precision k = MV n (A, A( n ~ m ^) to obtain a geometric resolution Ry of an 
equidimensional variety of dimension n — in. 

8. Obtain the geometric resolution R := A~ l Ry of V*(f). 
OUTPUT: The geometric resolution R of U*(f). 



In the sequel we will justify the correctness of the above procedure and estimate 
its complexity. 
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Since f is a generic sparse system of m equations in n variables, as a consequence 
of Lemma [fl if L%, . . . , L n _ m are generic linear forms, V*(f) D V(Li, . . . , L n - m ) is 
either the empty set (when V*(f) is the empty set) or a finite set consisting of 
deg(V*(f)) points (when V*(f) is not the empty set). Step CD decides whether V*(f)D 
V(Li, . . . , L n _ m ) is empty or not, since the mixed volume of the supports of these 
polynomials is the sum of the volumes of the mixed cells in a fine mixed subdivision 
[13]. If it is not empty, this finite set of points can be regarded as a generic fiber 
of a generic linear surjective projection and therefore, it enables us to recover the 
variety V*(f) by deformation techniques. 

Thus, the idea of the algorithm is to choose n — m linear forms at random, 
then compute a geometric resolution of the set V*(f) fl V(Li, . . . , L n _ m ) and finally, 
apply a Newton-Hensel lifting to the finite set obtained in order to get a geometric 
resolution of V*(f). 

Step [2] deals with the random choice of the entries of a matrix and a vector. 
We will not consider the cost of this random choice in the overall complexity of the 
procedure. In Step [3J the sparse encoding of n — m linear forms constructed from 
the previous data is obtained. 

The idea of Step 0] is to obtain a geometric resolution of V*(f)nV(Li, . . . , L n _ m ). 
In order to do this, the algorithm computes the isolated common zeroes in (C*) ra of 
the generic system f , L 1; . . . , L n „ m supported on A, A^ n ~ m \ Note that, if L 1; . . . , L n _ Ti 
are generic, this set of points meets only the irreducible components of V*(f), that 
is, it contains no point in the irreducible components of V(f) with vanishing coordi- 
nates. By applying the result in [TTf Proposition 5.13], it follows that the complexity 
of this step is O ((n 3 (N + (n - m)n) \ogd + n 1+n )M(D)M(m)(M(D) + M(E))), 
where 

. JV := £ 4Af, 
l<i<n 

• d := maxi< J < n {deg(/ j )}; 

• D := MV n (A, A (n ~ m) ); 

• 9Jt := max{||/i||}, where the maximum ranges over all primitive normal vectors 
to the mixed cells in the fine mixed subdivision of (A, A^ n_m ^) given by u; 

• E := MV n+ i(Ax{0},Ai(uji), . . . ,Am(uj m ),A(u) m+ x), . . . , A(w n )), where Aj{uj) 
(1 < j < m) and A(oji) (m + 1 < / < n) are, respectively, the supports of 
f, L u . . . , L n _ m lifted by uj. 

Now the algorithm lifts the geometric resolution of the zero-dimensional subset 
of V(f) obtained so far to a geometric resolution of the union of the irreducible 
components of this variety having a non-empty intersection with (C*) ra . In order to 
do this, we consider the change of variables given by Y = A.X and make this change 
of variables in the polynomials f and the geometric resolution already obtained 
(Steps [5] and |6]). The cost of computing A' 1 is 0(n n ) and the cost of computing 
the product of this matrix times a vector is 0(n 2 ). Then, we obtain an sip of 
length L = 0(n 2 + nlog(d)N) for the polynomials in g. Taking into account that 
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the degrees of the polynomials v i, . . . , v n are bounded by D, to write the geometric 
resolution in the new variables, we perform 0(n 2 D) operations. 

Note that (u>i(w), . . . , w n - m (u)) = b and (g(u), w n - m+l {u), . . . , w n (u)) is a ge- 
ometric resolution of the isolated points in V(g(b,Y n _ m+ i, . . . ,Y n )) corresponding 
to the isolated points in (C*) n of V(f, Li, . . . , L„_ m ). Now, the geometric reso- 
lution of V*(f) with respect to the linear projection given by Y consists of poly- 
nomials in Q[Yi, . . . ,Y n _ m ,u] having total degrees bounded D. Therefore, it suf- 
fices to compute the representatives of these polynomials in (Q[Yi, . . . , Y n _ m ]/ (Yi — 
61, ... , Y n - m — b n - m ) D+1 )[u\. To this end, in Step Owe apply successively the Global 
Newton Iterator from [TT] to the polynomials g, starting with the geometric reso- 
lution (q(u), w n - m+ i(u), . . . , w n {u)) obtained in Step El which can be regarded as 
a representative in (Q[Y"i, . . . , Y n _ m ]/ (Y 1 - b%, . . . , Y n _ m - b n - m ))[u], up to the re- 
quired precision D = MV n (A, A*" - "" 1 )). Using [TTJ Lemma 2] and encoding the 
elements of Q[Y U Y n _ m }/ (Y t - b u . . . , Y n _ m - b n „ m ) k as (k + l)-tuples of sip's 
(one sip for each homogeneous component), the complexity of Step [7] is of order 
0(( y mL + m n )M(D)D 2 ). 

Finally, the algorithm changes variables back in order to obtain the desired 
geometric resolution of V*(f), which adds 0(n 2 D) to the complexity. 

Taking into account the previous complexity estimates, we have the following 
result: 

Proposition 2 Let f = . . . , f m ) be a system of m < n generic polynomials in 
Q[Xi, . . . ,X n ] supported on A = (Ai,...,A m ). GenericToricSolve is a proba- 
bilistic algorithm that computes a geometric resolution of the affine variety V*(f) 
consisting of the union of all the irreducible components of V(f) that have non- 
empty intersection with (C*) n . Using the previous notation, the complexity of this 
algorithm is of order 

O (n 3 (N + (n- m)n) \og(d)M(D)(M(Wt)(M(D) + M(E)) + D 2 )) . 

3.2 Affine components 
3.2.1 Theoretic results 

This section is devoted to showing a combinatorial description of the equidimen- 
sional decomposition of the affine variety defined by a generic sparse polynomial 
system. More precisely, we prove combinatorial conditions on the supports of the 
polynomials that determine the existence of irreducible components of the different 
possible dimensions not intersecting (C*) n and give a combinatorial characteriza- 
tion of the set of linear subspaces where these components lie. This characterization 
enables us to give a combinatorial formula for the degree of these varieties. 

The following example shows that generic square sparse systems may define affine 
varieties containing positive dimensional components. It also shows that neither the 
mixed volume nor the stable mixed volume of the system are upper bounds for the 
degree of the affine variety defined: 
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Example 1 Consider a generic sparse system supported on A = (Ai,A2, A3) where 
A 1 = {(1, 1, 2), (1, 1, 1)}, A 2 = {(2, 0, 1), (1, 0, 1)} and A 3 = {(0, 2, 1), (0, 1, 1)}; 



aX^Xl + bX±X 2 X 3 = 
cX\X 3 + rfXiX 3 = 
eX\X 3 + fX 2 X 3 = 

The affine variety defined by the system has 5 irreducible components of degree 
I: {x 3 = 0}, {x x = 0,x 2 = ~i},{xi = ~,x 2 = 0}, {x ± = 0,x 2 = 0} and 
-!)}■ However, MV 3 (A U A 2 , A 3 ) = 1 and SM 3 {A U A 2 ,A 3 ) < MV^U 
{0},A 2 U{0},A 3 U{0})=4. 

Let f = (/1, . . . , f s ) be generic polynomials in the variables X = (Xi, . . . , X n ) 
supported on A = {A\, . . . , A s ) C (Z" ) s . 
For / C {1, . . . , n}, we define 

Ji = {3 e {1, • • • , s} I /jln ieJ {xi=o} ^ 0}, 

that is, the set of indices of the polynomials in f that do not vanish identically under 
the specialization Xi = for every i E I, and 

f/ = ((fj)i)je.Ji 

where, for a polynomial / G C[X 1? . . . , X n ], fj denotes the polynomial in C[(Xj)j^/] 
obtained from / by specializing Xi = for every i £ I. Namely, f/ is the set of 
polynomials obtained by specializing the variables indexed by I to in the poly- 
nomials in f and discarding the ones that vanish identically. We denote by Aj the 
support of (fj)i, by 717 : C n — > C n ~ #/ the projection iri(xi, . . . ,x n ) = (x^^i onto 
the coordinates not in / and by ipi : C n_ * 7 — > C n the map that inserts zeroes in the 
coordinates indexed by /. 

For an irreducible subvariety W of V(f) = {x e C n | = for every 1 < 

j < s}, let 

/ B r = {ie{l 1 ...,n}|^c{i i = 0}}. 

Lemma 3 Under the previous assumptions, let W be an irreducible component of 
V(f). Then dimW = n — - i^-hw Moreover tti w (W) is an irreducible compo- 
nent ofV(f Iw ) C C n -* Iw intersecting (C*) n ~# Iw . 

Proof: Without loss of generality, we may assume that I\y = {r + 1, . . . , n} and 
Jj w = {1, . . . , m} for some r > and m < n. Then, n := n Iw : C n — > C r is the pro- 
jection to the first r coordinates and f/ = (fi(xi, . . . , x r , 0), . . . , f m (xi, ■ ■ ■ ,x r , 0)), 
where is the origin of C n_r . 

Note that W = tt(W) x {0} and tt(W) C V(f/). If tt(W) C C C V(fj) for an 
irreducible component C of V(fj), it follows that W C C x {0} C V(f), with C x {0} 
an irreducible variety. Since W is an irreducible component of V(f), the equality 
W = C x {0} holds. This implies that n(W) = C is an irreducible component of 
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In addition, by the definition of Iw, we have that W D yO^^lxj ^ 0}J 7^ 0: 

otherwise, W C \Sj=\{ x j = 0}; which implies that 17 C {x,- = 0} for some 1 < j < r 
since 17 is an irreducible variety. Therefore, iv(W) H (C*) r 7^ 0. 

By Lemma [[J we conclude that ir(W) D (C*) r has dimension r — m and so, 
dim(17) = dim(vr(17)) = n - # J w - #J/ W . □ 

Now we will characterize the sets I C {1, . . . ,n} such that the irreducible com- 
ponents of 7(f 7 ) intersecting (C*) n- * / yield irreducible components of V(f). 

Proposition 4 Under the previous assumptions, let I C {1, . . . ,n}. Then 7(fj) fl 
(C*) n_#/ is not empty if and only if for every J C Jj, dim(^ , g j Aj) > #J and, m 
this case, V A (f/)n(C*)™~* 7 zs an equidimensional variety of dimension n — #1 — j^Ji. 
In addition, if W is an irreducible component of V(fj) fl (C*) n ~* / , we have that 
ipi(W) is an irreducible component of 7(f) fl [^^{xi 7^ 0} if and only if for every 

la i, #7+#j r >#/ + #j J . 

Proof: The first statement of the Proposition follows from Lemma [TJ 

Without loss of generality, assume that I — {r + 1, . . . ,n}. Let W be an irre- 
ducible component of 7(fj) fl (C*) r . 

Suppose that for a subset I C I the inequality #J + # Jj < #1 + #Jj holds. 
Assume I — {f + 1, . . . , n) for f > r. 

Note that if £ = (&,..., f r ) G V(fj), then (f, n _ r ) G 7(f) and, therefore, 
(£,0f- r ) G 7(fj). Thus, we may consider W x {0f- r } C 7(fj), which will be 
included in an irreducible component 17 of 7(fj). Taking into account that ij 
consists of j^Jj polynomials in n — #1 variables and applying Lemma [1] to 17 and 
f/, it follows that 

dim(17) >n-#I-#J J >n-#I-#J I = dim(W). 

We conclude that <fi{W) = 17 x {0 n __ r } C 17 x {0 n _ f } C 7(f) and, therefore, 
<£>j(W) is not an irreducible component of 7(f) fl {xj 7^ Vz ^ /}. 

Conversely, if <pi(W) = 17 x {0 n _ r } is not an irreducible component of 7(f), 
there is an irreducible component 17 of this variety such that W x {0 n _ r } C 17. The 
previous inclusion implies that C J. Assume 1^ = {f + 1, . . . , n) for f > r. Due 

to Lemma [3J 7Tj~(T7) is an irreducible component of 7(fj_) having a non-empty 
intersection with (C*) r . Therefore, 

n - #I W - = dim( 7 r / _(t7)) = dim(W) > dim 17 = n - #/ - #J 7 , 

and so, #1 + #J/ > #/^ + # Jj_ . □ 

As a consequence of Proposition HI we have that the irreducible components of 
7(f) C C n are contained in the linear subspaces f] ieI {xi = 0} associated to the 
subsets I C {1, . . . , n} in 

r = {/ C {1, . . . ,71} I VJ C Jj, dimiJ^Al) > #J; V/C /, #J 7 +#7> #J/+#/}. 
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Note that there may be sets 1% C I 2 C {1, . . . , n} both in T, as it can be seen in 
Example [U where the three sets {1}, {2} and {1,2} give irreducible components of 
the variety. 

If we write V*(fi) to denote the union of all the irreducible components of V(fj) 
having a non-empty intersection with (C*) n- * J , from Lemma [3] and Proposition 0J 
we deduce that, for every I G T, 

vitrei)) U w - 

W irred. comp. of V(f) 
such that Iytr=l 

We obtain the following characterization of the equidimensional decomposition of 
V(f) and, using Lemma [U a combinatorial expression for the degree of V(f): 

Theorem 5 Let f = (/i, •••,/«) be generic polynomials in n variables supported 
on A = (Ai, . . . , A s ) C (Z> ) s . For k — 0, . . . , n, let 14(f) be the equidimensional 
component of dimension k ofV(f). Then, using the previous notations: 

v k (f) = |J <Pi(V{fi)). 

lev, 
#l+#J z =n-k 

Moreover, deg(V(f)) = £ /er MV n (A*, A^-* 1 -*^). 



Example 2 Consider the following system of generic polynomials m QLXl, X 2 , X 3 , X4 



a,\X\X$ + a 2 X\X\ + a%X\X2X-t, + 04^2^3 — 
b x X x X 2 + b 2 X 1 Xi + 63X1X3X4 + 64X3X4 + 65X3X4 2 = 
c 1 X 1 X 2 X i + c 2 X 1 X 3 X 4 + c 3 X 2 X 3 + 04X2X3X4 = 
UiXx + d 2 Xl + rf 3 XiX 2 + d 4 X| + C/5X3X4 = 

Then, with the previous notation, Y = {0, {1, 2}, {1, 3}, {2, 3}, {2, 4}, {3, 4}} and 
then, 

. v (f) = ^(f)u{(o,o,S,-|),(-|,o,o,^),(-| + |f,-|,o,o)}.- 

— For I = we have 19 isolated solutions in (C*) 4 (this quantity is given 
by the mixed volume of the family of supports associated to the system) 

— For I = {1,2} we have: 



f {l,2} 



64X3X4 + 65X3XI 
1 ,/ V V 



d^X\ + C/5X3X4 



^(f{i, 2} ; 



64^5 64 

65^ ' h 



which gives the point (0,0, j^, -^). 
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For I = {2, 3} we have: 

I{2,3} = < , v , , v2 ; V (t{2, 3 }J " 



which gives the point (— 0, 0, |^jr)- 
For / = {3, 4} we nave: 

.61X1X2 + o 2 XiXf \ _ J7 rf i , Ms b i 

i{m ~ + d 2 xf + d 3 x,x 2 ; K(f{3 - 4}) -\r^ + M^'^ 



w/mc/z azves t/ie pozni + |J,-^,0,0). 

• V^(f) = {igC 4 |i 2 = 0,i 4 = 0,diX! + d 2 x\ + d A x\ = 0} 

— For I = {2, 4} we have: 

f{2,4} = {d 1 X 1 +d 2 Xl+d i Xl; V*(i{ 2A y) = {(x l7 x 3 ) \ diXi+d 2 xl+d 4 xl = 
which gives the curve {x 2 = 0, x 4 = 0, di£i + d 2 :Ei + ^4^3 = 0}. 

• V 2 (f) = {iGC 4 |ii=0,x 3 = 0}: 

— For I = {1,3} we have: 

f {1>3} = 0; 7*(f {1)3} ) = C 2 , 
which gives the plane {x\ = 0, 23 = 0}. 

Remark 6 In the case of a generic system f = (/1, . . . , f s ) in n variables such that 
A± = ■ ■ ■ = A s , we have that the sets I G T, I ^ 0, are all I C {1, . . . ,n} such that 
#Ji = and for every I C / ; > 0; and G T if and only «/dim(^4i) > s. 

Moreover, for every I G T, I 7^ 0, y>/(V"*(fj)) = {xj = Vz G J} and so, apart 
from the components intersecting (C*) n that correspond to I — (if$£T), the only 
irreducible components 0/ V(f) are linear subspaces ofC n . 

3.2.2 Algorithmic results 

According to Proposition HI the subsets / C {1, . . . , n} which yield irreducible com- 
ponents of the variety V(f) are the ones satisfying simultaneously 

1. VJC J 1, dim(^^j) > #J, 

2. V/' C /, #J r + #/' > #J/ + #/. 
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Now we present a possible (exponential) algorithm to obtain the sets 7 satisfy- 
ing condition [2] and the inequality #1 + # Jj < n, which is a necessary condition 
for the system f/ to have zeroes in (C*) n "# J , weaker but easier to check than con- 
dition [TJ Our algorithm to find all the affine components of V(f) (see Algorithm 
GenericAf f ineComps below) checks only among these sets whether condition [TJ is 
fulfilled or not by means of a mixed volume computation. 



Algorithm SpecialSets 

INPUT: A family of supports A= (Ax,...,A s ) C (Z%) s . 

1. P0 := min{n, s}. 

2. If = s, add (0, {1, . . . , s}) to an empty list T. 

3. For k — 1, . . . , n: 

For every 7 such that #7 = k: 

(a) Pj := mm{n, {Pj'}j/ c j,#i'=*-i, fc + #<M- 

(b) If Pi = k + #Jj, add (7, J/) to the list f . 

OUTPUT: The list T of all pairs of subsets (7, J/), J C {1, . . . , n} such that #7 + 
# J/ < n and for every 7 C 7, #7 + #7 7 > #7 + # J/. 



First, let us prove the correctness of this algorithm: 

Lemma 7 For every I C {1, . . . , n} ; Pj = min{n, #7 + #7 7 }. 

7c/ 

Proof: By induction on #7. 

For #7 = 0, since #7g = s, we have that P0 = min{n, #0 + #70}. 

Assuming the identity holds for every subset of cardinality k — 1, let 7 C 
{1, ...,n} with #7 = fc. Consider a proper subset 7 C 7. Then, there exists 
7' C 7 with #7' = A; — 1 such that 7 C I'. By the inductive assumption, 

P P = mm{n, #7 + #7 7 } 
id' 

and so, P// < #7 + #7 7{) . On the other hand, by the definition of Pj, we have that 
Pi < P r . Thus, Pi < #7 + #7 7q . □ 

Now we estimate the complexity of the algorithm. Let N = ^2j = i#A.j. For 
each 7 of cardinality k > 1, it takes fcA 7 " + #7/ + 1 operations to compute k + # Jj. 
Taking the minimum among k + 2 numbers takes /c + 1 comparisons. Thus, the 
complexity of Step [3a] is k(N + 1) + #7/ + 2. In Step [30, we add one comparison. 
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As we have to do this for each subset of {1, . . . , n} of cardinality k > 1 and for 
the empty set, the complexity is bounded by 2 + J2k=i Gt)(M-^ + 1) + s + 3) = 
2 + n2 n ~ 1 (N + 1) + (2 n - l)(s + 3) which is of order 0(nA2 n ). 

Unfortunately, the exponential complexity of the algorithm cannot be avoided, 
as the following examples show: 

Example 3 Let L 1; . . . , L n G Q[Xi, • • • , A n ] 6e generic affine linear forms. Con- 
sider the set of generic polynomials f\ = X\.L\, . . . , f n = X n .L n . In this case, if 
A = (Ai, . . . , A n ) is their family of supports, F — {I | I C {1, . . . , n}} and therefore 

#r = 2 n . 

One may think the exponentiality of the cardinal of the set T in this example arises 
from the fact that the variables are factors of the polynomials. However, this is 
not always the case as the following example shows. This example also shows how 
subroutine SpecialSets is useful to discard subsets which do not lead to affine 
components: in this case, the only element in T which does not correspond to a set 
in T is (0,{l,...,2n». 

Example 4 Consider generic polynomials f ± , . . . , f 2n G Q[-X"i, • • • , X 2n \, 

fj = ^ Cj k X 2 k~iX 2 k, j — 1, • • • , 2n, 

l<k<n 

all supported on the set A = {e 2 k-i + e 2 k \ 1 < k < n} where denotes the ith 
vector of the canonical basis ofM. 2n . Then, it is easy to see that, for any subset 
S C {1, . . . , n}, the set I s = {2k - 1 | k G S} U {2k \ k G {1, . . . , n} \ S} is in T, 
and that the sets Is are the only ones (together with the empty set) obtained as first 
coordinate of elements of F by the previous algorithm. Therefore, in this case, we 
have that the number of elements of the list F is 2 n + 1. 

In the following example, the usefulness of subroutine SpecialSets is more 
evident: 

Example 5 Let A = (Ai, . . . ,A n ) be a family of finite sets of (Z> ) n such that, 
for every 1 < j < n, there exists a nonnegative integer dj. such that dj v ei G Aj 
for every 1 < i < n. Then the output of subroutine SpecialSets in this case is 
r = {(0,{l,...,n})}. 

The following algorithm computes a family of geometric resolutions describing 
the affine variety defined by a generic system f with given supports A. 



Algorithm GenericAf f ineSolve 

INPUT: A sparse representation of the generic system f = (/i, . . . , f s ) of polynomials 
inQ[AV--,*n]. 
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1. Apply algorithm SpecialSets to the family of the supports A = (Ai, . . . , A s ) 
of f to obtain the list F of pairs of sets (I, Jj) with / C {1, . . . , n} such that 
#/ + #Ji < n and V/ C /, #7+ # J 7 > #7 + # J/. 

2. For every (J, J/) G T: 

(a) For j G J/, compute := {717(a) \ a E Aj such that = Vi G /}, 
and obtain the sparse representation of the system f/ supported on A 1 = 

(b) Apply algorithm GenericToricSolve to the sparse system f/ to obtain 
a geometric resolution R 1 (possibly empty) of the affine components of 
the set of solutions of f/ intersecting the torus (C*) n ~ #/ . 

(c) If R 1 ^ 0: 

i. Obtain the geometric resolution ipj^R 1 ) of the union of all irreducible 
components W of V(f) such that Iyy — I by adding zeroes to R 1 in 
the coordinates indexed by /. 

ii. If n - (#1 + = fc, add v?/(^ 7 ) to the list V fc . 

OUTPUT: A family of lists Vfc, < A; < n — 1, of geometric resolutions, each list 
either empty or describing the equidimensional component of V(f ) of dimension k. 



The correctness of this algorithm is straightforward from Proposition H] and The- 
orem |5j 

When applying algorithm GenericToricSolve in Step l2bl again, we do not 
take into account the cost of a pre-processing obtaining a fine mixed subdivision, 
which in particular decides whether the set I satisfies dim(^ jgJ ^4j) > #J for every 

J C J i and, therefore, it discards the sets I G T \ Y. For this reason, we will only 
consider the complexity of the computations for the sets I G V. This complexity can 
be estimated from the complexities of the subroutines applied at the intermediate 
steps (see Proposition |2]) and is of order 

0(^(n-#/) 3 (iV / +(n-#/-#J / )(n-#/)) \og(di)M(Dj) (M(Wlj) (M(Dj) + M(Ej)) + D 
/er 

where, for every / G T, Ni, dj, Dj, DJti and Ej are the parameters defined in the 
complexity of Algorithm GenericToricSolve, associated to the system fj. 

In order to estimate the overall complexity of the algorithm, note that Nj < 
N := J2j=i#Aj and d I < d := max 1 < J < s {deg(/ 7 )} for every I G Y. In addition, 
V := ^2i er Dj = degU(f). Note that, if w max is the maximum value of the lifting 
functions applied to the supports A 1 , for every I G T we have 

Ej < MU n „ #/+1 (A x {0}, [A] x {0,uw}) ieJl , (A x {0, a; max }) (n - #/ - #J ^) < 

< ((«-#/-# J;)MK- # ;(/, A (n-#/-#J,) )+ ^ MK-#/(C4^, A^* 7 -*^ 1 ))) . 

leJi 
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Then, if £ max := max /er {(n - #/ - # J^MV^A 1 , ^-#1-*^)) + 
E^MK-^a^^A^-* 7 -*'^ 1 ))} and 2K max := max^ri^}, taking into 
account the complexity of Step [TJ we have: 

Theorem 8 Let f = (fx, . . . , f s ) be a system of generic polynomials in Q[X\, . . . , X n ] 
supported on A = (A\, . . . ,A S )- GenericAffineSolve is a probabilistic algorithm 
that computes a family of lists Vk, < k < n — 1, of geometric resolutions, each list 
either empty or describing the equidimensional component ofV{f) of dimension k. 
Using the previous notation, the complexity of this algorithm is of order 

O (n2 n N + n 3 (N + n 2 ) \og{d)M{V) (M(9JT max ) (M(V) + M(u; max £ max )) + V 2 )) . 

4 Arbitrary sparse systems 
4.1 An upper bound for the degree 

The aim of this section is to show a bound for the degree of the affine variety defined 
by a square system of sparse polynomials which takes into account its sparsity. 

The following example shows that the degree of an affine variety defined by a 
generic sparse system with given supports is not an upper bound for the degree of 
the variety defined by a particular system with the same supports. One may think 
the problem arises from the presence of irreducible components not intersecting the 
torus either for the generic or the particular systems. However, this is not the case: 

Example 6 Consider the following system: 

X X X 2 -X 1 -X 2 + 1 = (X 1 - 1)(X 2 - 1) = 
Xl X 3 -X 1 -X 3 + 1 = (X 1 -1){X 3 -1) = 
X 2 X 3 - X 2 - X 3 + 1 = (X 2 - 1) (X 3 - 1) = 

The variety defined by the system consists of 3 lines, each having a non-empty in- 
tersection with (C*) 3 . However, if A = (Ai, A 2 , A 3 ) is the family of the supports of 
the polynomials, the degree of the variety defined by a generic system with the same 
supports is MV 3 (Ai, A 2 , A 3 ) = 2. 

Our bound for the degree, which is stated in the following theorem, is the mixed 
volume of a family of sets obtained by enlarging the supports of the polynomials 
involved. 

Theorem 9 Let f = (f\, . . . , f n ) be n polynomials in C[Xi, . . . , X n ] supported on 
A = (Ai, . . . , An) and let V(f) = {x e C n | fi(x) = for every 1 < i < n}. Let 
A = {0, ei, . . . , e n } where e« the ith vector of the canonical basis ofM. n . Then 

deg(V(f)) < MV n (Ai U A, . . . , A n U A). 
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Before proving the theorem, we will fix some notation and definitions. 

Let Tj = i^{Aj U A) — 1 for 1 < j < n. Consider the morphism of varieties 
<p : C" ->• FxF - x-xf" defined by <p(x) = ((1 : x), (x a ) a€AlUA , . . . , (x a ) aeAiUA ). 
Let X = ip(C n ). For 1 < j < n, we denote by Lj the linear form in P rj given by 
the coefficients of the polinomial fj, that is to say, if fj = ^2 aeAj Cj A X a , then 

For each integer k (0 < k < n) and each subset S C {1, . . . , k} we define the 
variety X kt $ recursively in the following way: 

1- <^o,0 = X . 

2. Provided X kt $ is defined for every S C {l,...,k}, we define with 
T C {1, . . . , k + 1} as follows: 

• If k + 1 ^ T, Xfc + i 5 T is the union of the irreducible components of X^x 
included in {L k+ i = 0}. 

• If k + 1 G T, Xk + i t T is the intersection of {L k+ i = 0} with the union of 
the irreducible components of -X^yyo+i) not included in {L k+1 = 0}. 

Note that, from this definition, each X k ,s is an equidimensional variety of dimen- 
sion n — #S. Moreover, if ir : P n x P n x ■ ■ • x P r ™ — y P n is the projection onto the 
first factor, it is easy to see inductively that, for every 1 < k < n, 

|J n(X k , s ) = Hfi,...,/*) C P n . 
sc{i,...,fe} 

For an equidimensional subvariety W of X, we define its multidegrees deg^ r 0kln _^(W) , 
for r, A; G Z> such that n — k + r = dim(W), as 

{r n ^ 

#(wnf|{4w = o}n f| {4 = 0}) 

where the maximum is taken over all (£q,i, ■ ■ ■ , io,r, 4h-i> ■ ■ ■ , ^n) such that each £ j 
is a linear form in P n and each lj is a linear form in P rj and the intersection is 
finite. Note that the 1 subscript indicates how many projective spaces are cut by a 
linear form and the subscript shows how many remain uncut. As in the case of the 
standard degree of affine or projective varieties (see p2]), the maximum is attained 
generically. 

In particular, it is clear that deg( 0i00)ln ) 

Lemma 10 Under the previous assumptions, definitions and notations, 

de S(k-#s,o k ,i n - k )( x k,s) > deg (fc+1 _ #S)0fe+1)ln _ fc _ l) (A' A+1|g )+deg (A ._ #5)0) , +iiln _ fc _ i) (^ +liSu{fe+1} ) 

Proof: As the variety X k> s = X k+ i^X k ^s, where X kj s is the union of the irreducible 
components of X k ,s not contained in {L k+ i = 0}, using genericity in the definition 
of multidegrees, we have that 

deg( fc _ #S) o fc ,i n _ fc )(^ 1 s) = deg (k _ #SfikAn ^ k) {X k+hS ) + deg {k _# Sfikjln _ k) {X kiS ). 
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Note that the points of the varieties we are considering are of the form 

(x, (x : (x a ) aeAl \ A ), ...,(x: (x a ) aeAn \ A )) . 

Because of this, if is a linear form in P n and ik+i is exactly the same linear 
form involving only the first coordinates in P rfc +!, we have that Xk+i,S H {io = 0} = 
Xk+l,S H {ik+i = 0} and therefore, 

de g(k-#s,o k ,i n _ k )( x k+i,s) > deg (fe+1 _ #Sj0fc+liln _ fc _ l) (A' fc+liS ) 

and 

deg (k _ #St0k)ln _ k) (X k>s ) > deg (fc _ #s>0jH . liln _ k _ l) [X ktS n {L fc+1 = 0}) = 

= de g(fc-#5,0 fe+1 ,l n -fe-i)( ; ^+ 1 .SU{fe+l})- 

□ 

Now, we can prove Theorem [9j 

Proof: Consider the projection ir : P n x P ri x • ■ • x P r ™ — > P n onto the first factor. 

As ' ' 

|J n(X n>s ) = V(f) C F\ 

Sc{l,...,n} 

we have that 

deg(V(f)) = deg ( |J 7r(A^)) < ^ deg7r(* n , s ) = 

Sc{l,...,n} Sc{l,...,n} 

= de g(n-#5,0„,l )( ;t 'r l ,5) 
5c{l,...,n} 

(this last equality is nothing but our definition of multidegree) . 

Applying inductively Lemma fTUl we get that, for each < k < n, 

de g(k-#s,o k ,i n . k )( x k,s) < deg (0j0o ,i„) X o,$ 

Sc{l,...,fc} 

and therefore we obtain that 

deg (n _ #Sj0n)lo) (A; )S ) < deg (o>OO)ln) (Ar o>0 ) = MV n {A x U A, . . . , An U A). 

Sc{l,...,n} 

□ 

In the following examples, we show that this bound may be attained: 
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Example 7 Let S C {1, . . . ,n} and let fx, . . . , f n G Q[Xi, . . . , X n ] be polynomials 
of degree d in the variables (Xiji^s with no zero coefficients (that is to say, the 
monomials appearing in fi, . . . , f n are those of degree less or equal to d in the vari- 
ables (Xi) ie s). If A is their common support, it is evident that MV n (A^) = and 
that SM n (A (n) ) = 0. However, MV n ((A U A) (n) ) = d# s and this degree can be at- 
tained for special choices of the polynomials: If fx, . . . , f#$ £ Q[PQ)ies] ^ s a family 
with d# s isolated solutions in C* s , take generic linear combinations f#s+i, ■ ■ ■ , fn °f 
fi, ■ ■ ■ , f#S- Then, the family of polynomials fx, ■ ■ ■ , f#s, f#s+i, ■ ■ ■ , fn e Q[X U ...,X n ] 
defines a variety formed by d# s affine linear spaces of dimension n — #S. 

Our bound can be also attained for generic systems: 

Example 8 Consider the family of generic polynomials fx, ■ ■ ■ , fn defined in Exam- 
ple\M and let their supports A = {Ax, ■ ■ ■ , An). Then, the affine variety they define 
consists of 2 n points, and therefore, its degree is 2 n = MV n (Ax U A, . . . , A n U A). 

4.2 An algorithm in the non-generic case 

In the sequel we will describe an algorithm that, given an arbitrary square system 
of sparse polynomials, provides a set of finite points in each irreducible component 
of the affine variety the system defines. The complexity of this algorithm depends 
on the bound for the degree of the variety obtained in the previous section. 

Let f = (fx, ■ ■ ■ , fn) be n polynomials in Q[Xx,...,X n ] supported on A = 
(Ax, ■ ■ ■ , A n ) and, for a fixed k, 1 < k < n — 1, let Lx, ■ ■ ■ , Lk be generic affine 
linear forms in Q[Xx, ■ ■ ■ , X n ]. Then, if Vfc(f) is the equidimensional component 
of dimension k of V(f), the isolated common zeroes of fx, ■ ■ ■ , f n ,L\, ■ ■ ■ ,Lk are 
deg(Vfc(f)) points in 14(f). The idea is to represent each equidimensional compo- 
nent by means of the corresponding set of points (c.f. the notion of witness point 
set in [30]). 

Example 9 Consider the following polynomial system: 

' XfX 2 X 3 - XxX 2 Xl - X\ + X\ = (XxX 2 X 3 - l)(Xx - X 3 )(Xx + X 3 ) 

„ _ XfX 2 X 3 — X\X 2 X 3 — XxX 2 X 3 + XxX 2 X 3 — XxX 2 + Xx + X 3 X 2 — X 3 = 
| = (XxX 2 X 3 -l)(Xx-X 3 )(X 2 -l) 

, XxX 2 X 3 — XxX 2 X 3 — X 2 + X 3 = (XxX 2 X 3 — 1)(X 2 — X 3 ) 
The equidimensional components ofV(f) are 

Vb(f) = {(-1, 1, 1)}, Vx(f) = {xx-x 3 = 0, x\-x 3 = 0}, V 2 (i) = {xxx 2 x 3 -l = 0}. 
Taking Lx = Xx — X 2 and L 2 = 6X 2 — X 3 + 7, we have: 

• the set of isolated points of V(f) fl {xx — x 2 = 0} is {(1, 1, 1), (0, 0, 0)}, which 
is a set with 2 = deg(Vi(f)) points in Vx(f). 

• V(f) H {xx - x 2 = 0, 6x 2 - x 3 + 7 = 0} = {(-1, -1, 1), (-i, -\, 4), (|, \, 9)}, 
which is a set with 3 = deg(V2(f)) points in ^(f). 
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For a fixed k, < k < n — 1, in order to compute the isolated common zeroes of 
fx, . . . , f n , Li, . . . , Lk, by taking n generic linear combinations of these polynomials, 
we obtain a system of n polynomials in n variables having these points among its 
isolated zeroes (see [12J). Note that, in order to achieve this, it suffices to take linear 
combinations of the form 

k 

f t (X) + Y,hMX), z = l,...,n 

i=i 

for generic bij (1 < i < n, 1 < j < k). 

Procedure PointsInEquidComps described below computes a family of n geo- 
metric resolutions RP°\ for < k < n — 1, encoding a finite set of points and such 
that RS k ' represents at least deg(J¥) points in each irreducible component W of 
dimension k of V(f). 

The intermediate subroutine CleanGR takes as input a geometric resolution 
(q(u), V\ («),..., v n (u)) C (Q[u]) n+1 of a finite set of points V C C n and a list of 
polynomials f = (fx, ... , f n ) G Q[Xi, . . . , X n ], and computes a geometric resolution 
(Q(u),V 1 (u),...,V n (u))oiVnV(f): 

Q(u) = gcd (q(u), fx {v x (u),...,v n (u)),..., /„(^i («),..., f n («))) 
Vi(ti) = fj(-u) mod Q(m), i = l,...,n. 

Let A,- be the support of fj, d an upper bound for deg(/j), j = l,...,n, and 
D = degq. First, the subroutine computes sip's of length 0(nD log D) for the 
polynomials Vi, i = l,...,n. The gcd Q{u) is computed successively as follows: 
For j = l,...,n, the subroutine computes an sip of length Lj = 0(n\ogd#Aj) 
for the polynomial fj and, by multipoint evaluation and interpolation, the dense 
representation of Fj(u) = fj(v i(w), . . . , v n {u)) within complexity 0(M(dD)(Lj + 
nDlogD)); then, it obtains Qj(u) := gcd(Qj-x,Fj(u)) within 0(M(dD)) addi- 
tional operations. Finally, the polynomials Vi(u) for i = l,...,n are obtained 
within complexity 0(nM(D)). The overall complexity of the procedure is of or- 
der 0(M(dD) (n log d YTj=i #Aj + n 2 D\ogD)). 



Algorithm PointsInEquidComps 

INPUT: A sparse representation of a system f = {fx, ■ ■ ■ , fn) of polynomials in 
Q[Xi, . . . , X n ] supported on A = {Ax, ■ ■ ■ , A n ), a lifting function uj = {cox, ■ ■ ■ , oj n ) 
for Aa = {Ax U A, . . . , A n U A) and the mixed cells in the induced subdivision of 

1. Choose randomly coefficients for a polynomial system g = {gx,...,g n ) sup- 
ported on Aa- 

2. Apply the algorithm in [TTt Section 5] to g to obtain a geometric resolution 
R s of its zeroes in C n . 
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3. Choose randomly n — 1 affine linear forms L±, . . . , L n _\ in the variables X = 
(Xi, . . . ,X n ) and n(n - 1) integer numbers (%)i<i<n,i<j<n-i- 

4. For k = 0, . . . ,n - 1: 

(a) Obtain the sparse representation of the polynomials h[ k \x) = fi(X) + 
X^ =1 bijLj(X) for 1 < z < n. 

(b) Apply the algorithm in pU Section 6] to h (fc) = (/if } , . . . , h { n ] ) to obtain 
from R s a geometric resolution of a finite set Vk which contains the 
isolated common zeroes of h (fc) in C n . 

(c) Apply subroutine CleanGR to the previous geometric resolution and f to 
obtain a geometric resolution of Vk nV(f). 

OUTPUT: The n geometric resolutions for < k < n — 1. 



In the sequel we will estimate the complexity of this procedure. 

Steps [T] and [3] are fulfilled by taking a random choice of numbers. We will not 
consider the cost of this random choice in the overall complexity. The complexity 
of Step [2] is 

0((n 3 N A log(d) + n 1+n )M(D A )M(M A )(M(D A ) + M(E A ))) 

where 

• AT a :=E; =1 #(A-UA); 

• d := maxi^^deg/j-}; 

• D A := MV n (Ai U A, . . . , A U A); 

• 9JTa := max{||/i||} where the maximum ranges over all primitive normal vec- 
tors to the mixed cells in the fine mixed subdivision of A A induced by u; 

• E A := MV n+1 {A x {0}, (A U A)^), ...,{A n U A)(w n )) where (A, U A)(w i ) 
for every 1 < j < n is the set A,- U A lifted by u>. 

In Step Ebl we compute a finite set which includes the affine isolated zeroes of 
the system h^. By applying the result in [TTJ Proposition 6.1], we have that the 
complexity of this step is bounded by 0((n 2 N A logd + n 1+n ) M \D A ) M (E' A )) where 
E' A := MU n+ i({0} x A, {0, 1} x (A% U A), . . . , {0, 1} x (An U A)). Finally, the 
complexity of Step He] is of order 0(M(dD A )(n log ^EJ=i #A' + n 2 D A \og D A )). 

Note that the parameters E A and E' A in the previous complexities can be bounded 
as follows: 

n 

E' A = ^2 M K(A, Ail) A,..., AfuA, . . . , An U A) < nD A 

3=1 

and, if w max : = max -j,a{wj(a) | 1 < j < n, a G Aj U A}, 

E < MV n+ i(A x {0}, (Ai U A) x {0,w max }, . . . , (A n U A) x {0, w max }) < cu mBX nD A . 
Taking into account these bounds, we have: 
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Theorem 11 Let f = (/i, . . . , f n ) be n polynomials in Q[Xi, . . . , X n ] supported on 
A = (Ai, . . . , An). Point slnEquidComps is a probabilistic algorithm which com- 
putes a family of n geometric resolutions (R^°\ . . . , R^ n ~^) such that, for every 
< k < n — 1 ; R( k > represents a finite set containing deg Vfc(f) points in the equidi- 
mensional component Vfc(f) of dimension k ofV(f). Using the previous notation, 
the complexity of the algorithm is of order 

0(to miiX n 4 N A \og(d)M(dD A )M(D A )M(Tl A )). 

Example 10 Consider the polynomial system introduced in Example 0, given by 
the polynomials 

fx = X1X2 — X\ — X2 + 1 
f = lf 2 = -X 1 -X 3 + l 

% /3 — x 2 x 3 — x 2 — x 3 + 1 

supported on A 1 = {(1, 1, 0), (1, 0, 0), (0, 1, 0), (0, 0, 0)}, A 2 = {(1, 0, 1), (1, 0, 0), (0, 0, 1), (0, 0, 0)} 
and A 3 = {(0, 1, 1), (0, 1, 0), (0, 0, 1), (0, 0, 0)} respectively. 

Algorithm PointsInEquidCornps first chooses (at random) a system supported 
on (Ai UA,i 2 U A, i 3 U A), for example: 

2X 1 X 2 -2X 1 + X 2 -X 3 + 1 
X X X 3 -X 1 + 2X 2 + 2X3 + 2 
X 2 X 3 + X\ — 2X2 + X 3 — 1 



and computs a geometric resolution R g of its isolated common roots in C 3 : 



Re 



155 



u 5 - fn 4 - 17u 3 + 80u 2 -2u- , 







x x = + > 3 + IV -Mu- 137 



x 2 = -^M 



100 

200 



200 



200 



10 
7 



80 400 ^ 400 ~ 80 



Xo = -± U 4 + — "U 3 + V - ^-U - i3 
- 3 10 ^ 20 ^ 4 20 4 

Then, in Step® the algorithm takes 2 linear forms: 



L l = X t +X 2 + 2X 3 
L 2 = X\ + 2X2 

In Step\Jj^ for k = 0, 1, 2 ; i/ie isolated roots of the system h.^ obtained by adding to 
f generic linear combinations of Li, i — 0, . . . , k, are computed: 



f 1 (X) + L 1 (X) 

h (°)=f, h« = <! / 2 (I)-i!(I) 

/ 3 (X) + 2L 1 (X) 

/ 1 (X) + L 1 (X) + L 2 (X) 
h( 2 ) = <! / 2 (X)-L 1 (X) + 2L 2 (X) 
/ 3 (X) + 2L 1 (X) + L 2 (X) 



: X X X 2 + 2X3+1 

X X X 3 - 2Xi - 3X 3 - X 2 + 1 
: X 2 X 3 + X 2 + 3X 3 + 2X 1 + 1 

XxX 2 + X t + 2X 2 + 2X3 + I 
XxX 3 + 3X 2 -3X3 + I 
X 2 X 3 + 3X 2 + 3X 3 + 3X X + 1 
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In order to do this, the algorithm deforms the geometric resolution R g to geometric 
resolutions R h (k) of the sets of isolated roots o/h^: 



' u 3 



R h (o) = < 



7m 2 + 2u + 40 = 
X 1 = l 



V — 1 , ,2 i 9 3 
A 2 — -jjU + Yi M 7 



X, 



14 



-Rht 1 ) 



m 5 - f w 4 - 13m 3 + 68m 2 - 64m = 

57 
X 

_ _269 



y- _ _5Y_ 7 ,4 _ 369 7 .3 _ 953„.2 , 647„, _ o 
^ 1 200 " 400 " 200 " ^ 50 " d 



-2 

x 3 



,A I 1573„.3 , 1567 7( 2 _ 2599 , i 

600 " ^ 1200 " ^ 200 " 150 " ^ 1 



367 4 
600 



2039 3 
1200 



2181 2 , 3407 , 1 
200 " ^ 150 " ^ X 



i? h ( 2 



M° + f M 4 - ^U 6 + + - ^ = 



1549 „, 3 



538„,2 



679, 



769 



9 



x x = 


101214 4 
1803049 


x 2 = 


58338 4 , 
9015245 ~ l ~ 




389394 4 , 
9015245 u ^ 



9 

2537721 
1803049 



U 



18 

, 15987545 2 
i ~ 1803049 

11347628 2 i 
9015245 ~ l ~ 



277533 3 
1803049 



1982655 3 _ 57242469 2 _ 

I" nm r:o/i c; ^ 



I 3986650 
^ 1803049 
5343077 



U 



1803049 ' 



9015245 



9015245 

21604159 
9015245 



U + 

u4 



_ 7719426 
1803049 
8036523 
9015245 
22524084 
9015245 



Finally, subroutine CleanGR removes spurious factors from R h (k) to obtain a geo- 
metric resolution R^ of a finite set containing a set of representative points of the 
equidimensional component of dimension k 



RW 



u 3 - 7u 2 + 2u + 40 
X l = l 







X, = Tjlt 2 



14 



—U — — 

14" 7 



I u 2 _ 2 _ 1 
7 U 7 U 7 



m 3 + 2m 2 - 8w = 
Xi = \u 2 + U - 3 
X 2 = -\u 2 + \u+ 1 
X 3 = - V - |u + 1 



it>( 2 ) 



Since u 3 — 7m 2 + 1u -\- 40 = (it + 2) (u- 4) (it -5) and u 3 + 2u 2 - 8u = (it + 4)u(u-2), 
substituting their roots in R^ and R^ respectively, we get the following sets of 
points in V(f): 

W = {(1, -2, 1), (1, 1, 1), (1, 1, 2)} Wi = {(-3, 1, 1), (1, 1, -1), (1, -3, 1)}. 



Note that Wi contains exactly one representative point for each of the 3 lines in 
V(f). Moreover, the fact that R^ = implies that V(f) does noi /lave irreducible 
components of dimension 2. However, although there are no isolated points in V(f), 
Wq C V(f) noi empty. 
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